CMB multipole measurements in the presence of foregrounds 
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Most analysis of Cosmic Microwave Background spherical harmonic coefficients aim has focused 



on estimating the power spectrum Ci 



) rather than the coefficients themselves. We present 



a minimum- variance method for measuring aim given anisotropic noise, incomplete sky coverage 
and foreground contamination, and apply it to the WMAP data. Our method is shown to consti- 
tute lossless data compression in the sense that the widely used quadratic estimators of the power 
spectrum Ci can be computed directly from our a^m-estimators. As the Galactic cut is increased, 
the error bars Aaf,„ on low multipoles go from being dominated by foregrounds to being dominated 
by leakage from other multipoles, with the intervening minimum defining the optimal cut. Applying 
our method to the WMAP quadrupole and octopole as an illustration, we investigate the robustness 
of the previously reported "axis of evil" alignment to Galactic cut and foreground contamination. 



I. INTRODUCTION 

Cosmology has been revolutionized by the advent of 
precision maps of the Cosmic Microwave Background 
(CMB), allowing accurate tests of the Cosmological Stan- 
dard Model (CSM) and measurements of its free param- 
eters {e.g., [1-5]). Although the analysis of CMB maps 
has traditionally focused on the power spectrum, various 
surprises discovered in the data have triggered a grow- 
ing interest in extracting additional information pertain- 
ing to possible non-Gaussianity. Specifically, the all-sky 
CMB temperature fluctuation map 6T{r) is customarily 
expanded in spherical harmonics: 



(5T(f) =^a,™yft„(?). 



(1) 



These multipole coefficients aim are treated as stochastic 
variables, varying randomly when extracted from CMB 
maps in widely separated Hubble volumes. According 
to the CSM, the CMB fluctuations are for all practical 
purposes Gaussian and statistically isotropic, which im- 
plies that the a^m-coefficients are independent Gaussian 
random variables with zero mean, i.e., that the CMB 
contains no cosmological information whatsoever except 
for its power spectrum Ce = (la^mP)- 

Although some inflation models predict potentially ob- 
servable departures from Gaussianity on small angular 
scales {£ > 10^), no such non-Gaussianity has yet been 
found (see [6, 7] and references therein). However, a num- 
ber potential non-Gaussian anomalies have recently been 
discovered on large angular scales. The surprisingly low 
CMB quadrupole C2 has intrigued the cosmology com- 
munity ever since it was first observed by COBE/DMR 
[8], and was taken more seriously once the precision 
measurements of WMAP [9] showed that it could not 
be complete blamed on Galactic foreground contamina- 
tion. Although the low quadrupole does not involve non- 
Gaussianity and its statistical significance is debatable 
[1, 10, 11, 44], it is a fly in the ointment of the CSM, and 
closer investigation of this has revealed numerous hints 
of non-Gaussianity on large angular scales. 



The first reported hints of such non-Gaussianity in- 
volved the quadrupole and octopole. They were both 
found to be rather planar, i.e., with most of their hot and 
cold spots centered in a single plane in the sky, with their 
two preferred planes surprisingly closely aligned [12, 13]. 
In the CSM context, the reported alignment required a 1- 
in-66 fluke and the octopole planarity an independent 1- 
in-20 coincidence [13]. These features were confirmed by 
other groups and examined with other techniques, find- 
ing puzzling alignments all the way up to ^ = 5 that may 
require as much as a 1-in-lOOO fluke [13, 15-17, 20, 24- 
26, 28-30, 32, 49]. The preferred axis, dubbed the "axis 
of evil" by [25] , points towards Virgo and is intriguingly 
close to the ecliptic poles [16, 30, 49]. Additional hints 
of low-i non-Gaussianity have also turned up involving, 
e.g., asymmetries [14-42]. 

These purely observational results have triggered nu- 
merous papers on potential physical explanations. These 
include effects of Galactic foreground emission [43-47] 
and local structures [48-56] as well as theoretical expla- 
nations going beyond the CSM and involving compact 
cosmic topology [13, 57-69], modified inflation [70-83], 
or other new physics [84-101]. 

A substantial fraction of this work has involved the 
Ofoji-coefficients from equation (1) with £ < 5. Given this 
interest in the low-£ multipoles, it is timely and worth- 
while to measure them as accurately as possible, with 
quantified error bars, further improving on approaches 
such as [12, 18]. This is the purpose of the present pa- 
per. More accurate a£rn,-measurements can either build 
confidence in the CSM or teach us about interesting new 
physics. 

Because of Galactic foreground contamination and 
anisotropic detector noise noise, the best way to mea- 
sure the multipole coefficients aim is not to simply apply 
the relation aim = J Yem{^)*ST{r)dn to a CMB map 
(see Figure 1). Rather, the customary approach is to 
discard the most contaminated part of the sky and use a 
more elaborate linear weighting on the remainder of the 
map. In Section II, we study how to optimize this pro- 
cedure. We then apply our method to the latest WMAP 
data in Section HI, investigate the implications for the 



FIG. 1: Foreground contamination is clearly visible in the WMAP V-band map (top left). After taking a simulated WMAP CMB 
map (top right), adding simulated noise and foregrounds, applying the TOH cleaning method [12] and subtracting the CMB back out, 
foregrounds remain clearly visible in the residual map (bottom left). We apply our multipole measurement technique after masking out 
each of seven regions of decreasing cleanliness (bottom right) to optimize the tradeoff between residual foreground contamination and 
limited sky coverage (which causes leakage from unwanted multipoles coupling to our estimator). 



quadrupole-octopole alignment in Section IV and sum- 
marize our conclusions in Section V. 



II. BASIC METHOD 

A. The problem 

Given a CMB map with pixels i = 1, ...,N, let Xi de- 
note the observed temperature fluctuation ST in the pixel 
whose sky direction corresponds to the unit vector . For 
the case of WMAP, there are iV = 12 x 512^ = 3,145,728 
pixels distributed according to the HEALPix scheme^ 
[102, 103], but we will only use those outside of some 
appropriate Galaxy cut below, and at a lower resolu- 
tion. Suppose we wish to measure a particular set of 
multipole coefficients ae^m,, j = 1,...,M, labeled sim- 
ply aj for brevity^. Grouping the pixels into an iV- 



^ The HEALPix package is available from http://www.eso.org/ 
science/healpix/. 

^ For convenience, since CMB maps are real rather than complex- 
valued, we work with real-value spherical harmonics Ytm 
throughout this paper. These are obtained from the stan- 



dimcnsional vector x and the multipole coefficients into 
an M-dimensional vector a, we can rewrite equation (1) 
as 

X = Ya + n, (2) 

where the N x M spherical harmonic matrix Y is de- 
fined by Yij = Yg.jnji^i) and the "noise" vector n con- 
tains all contributions to the sky map except from the 
desired multipoles. In other words, n contains not only 
detector noise, but also genuine sky signal contributed by 
other multipole coefficients that are not included in the 
a-vector. We make the usual assumption that the noise 
has zero mean ((n) — 0) and define its covariance matrix 
C = (nn*). If the detector noise covariance matrix is N, 



dard spherical harmonics by replacing e"""^ by \/2 sin m(/>, 1, 
\/2cosm0 for m < 0, m = 0, m > 0, respectively. As detailed 
in Appendix A of [13], the corresponding real- valued coefficients 
a£j^ are related to the traditional complex-valued coefficients 
by a simple unitary transformation: for m < 0, m = and 
m > 0, equals \/2Ima£^, a^o ^'^d V2Reaf„, respectively. 
This means that for m < 0, m = and m > 0, the traditional 
complex coefficients equal (— l)'"(afm ~ ^'^e,-m) / V^, ciim 
and {ai,„ + ia^ _„)/V2, respectively. 
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then we can write this in the usual form 

C = N + S, Sij=Y,^^Pe{ri-rj)Ci, (3) 

e 

where S is the CMB contribution, Pi is a Legendrc poly- 
nomial and and the sum runs over those multipoles that 
are not included in our a- vector. We ignore the issue 
of foreground contamination for now, but will cover this 
important issue in detail below in Section III. 

We wish to find an estimator a of the true multipole 
vector a that is unbiased ((a) = a) and whose elements 
have as low variance as possible. 

B. The solution 

For our \ow-£ applications, the linear system given 
by equation (3) is normally greatly overdetermined with 
N ^ M, i.e., with many more pixels than desired multi- 
pole coefficients. Mathematically equivalent linear prob- 
lems frequently occur in the CMB literature in the con- 
texts of mapmaking and foreground removal, and the op- 
timal solution is well-known to be [104] 

a = Wx, W= [Y*C-^Y]-^Y*C-\ (4) 

with covariance matrix 

S^(aa*)-(a)(a)* = [Y*C-iY]-^ (5) 

To make this method useful in practice, we need to 
augment it with a prescription for how to optimize the 
Galaxy cut when faced with foreground contamination 
— we will do this in Section III. Before that, let us now 
provide some intuition for how the method works. 

First note that since the spherical harmonics are or- 
thogonal over the full sphere, complete sky coverage with 
uniform detector noise per pixel will make the covariance 
matrix S diagonal, reducing our method to simply ex- 
panding the map in spherical harmonics the usual way. 
The measurement of a given harmonic aem is then the 
exact sky signal plus a detector noise contribution. 

If part of the sky is missing due to foreground masking 
or lack of observation, or if the pixel weighting varies 
because of nonuniform detector noise, then the sphe- 
rical harmonics are no longer orthogonal, which is re- 
flected in S from equation (5) having off-diagonal ele- 
ments. However, the matrix algebra corrects for this cou- 
pling between the multipole coefficients. Equation (4) 
implies that WY = I so that the measurement error 
a — a = W(Ya -I- n) — a = Wn, independent of a. For 
example, if a contains the 16 multipole coefficients with 
^ < 3, then the measured quadrupole coefficient 020 will 
equal the true full-sky value 020 plus a noise contribu- 
tion. This "noise" will include detector noise and leak- 
age from sky multipoles with £ > 3, but there will be 
no contribution whatsoever from the monopole, dipole 
or octopole, nor from quadrupole coefficients a2m with 



TO 7^ 0. The method will generically be able to solve for 
these 16 unknowns as long as there are at least 16 pixels, 
but the error bars will clearly grow as the Galaxy cut is 
increased, since this makes it more difficult to disentan- 
gle the different multipoles. Numerically, we indeed find 
that the matrix [Y*C~'^Y] remains nonsingular for all 
the Galaxy cuts we consider, so no regularization tech- 
niques are needed. However, we will see that the error 
bars grow sharply as the cut increases and the matrix 
becomes less well-conditioned, particularly because mul- 
tipoles with large \m\ live predominantly in the masked- 
out regions near the Galactic plane. 



C. Relation to other methods 

The multiple estimation paper most closely related 
to this one is that of [18]. They consider two alterna- 
tive linear techniques that, cast in the notation of the 
present paper, simply use different W-matrices than the 
one given by equation (4) , perform a careful and detailed 
study of their statistical properties, and apply them to 
the WMAP data. The first technique they explore is 
Wiener filtering, defined by W = Y*S-i[YtC-iY]-i. 
Using numerical simulations, they confirm that this 
causes a systematic suppression of power with increas- 
ing £ (decreasing signal-to-noise) as expected. They then 
focus on the power equalization (PE) filter defined by 
a W-matrix based on Cholesky decomposing the sphe- 
rical harmonic coupling matrix Y*Y. The PE filter has 
the attractive property of eliminating all contributions to 
a given multipole estimator from lower multipoles, but 
it is neither unbiased (in the sense that WY 7^ I) nor 
minimum- variance. For the special case where we include 
all multipoles up to some £ in the a-vector, our estima- 
tor for the last a^m-coefficient will also be independent of 
the lower multipoles, merely with a lower variance than 
the PE estimator (since the estimator defined by equa- 
tion (4) by construction gives the smallest variance of 
any estimator with this property). 



D. Relation to power spectrum estimation 

The most common technique in the CMB commu- 
nity for measuring CMB power spectra uses so-called 
quadratic estimators [105, 106], largely because they 
have been shown to be information-theoretically optimal, 
giving the smallest possible error bars. Equations (4) 
and (5) imply that 

a*S-2a = x*C-iYY*C-ix. (6) 

When Y contains all spherical harmonics for a given 
£, the right hand side is precisely the quantity that 
a quadratic estimator of Ce extracts from the map x. 
Equation (6) therefore shows that we can calculate opti- 
mal estimators of the power spectrum from our measured 
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multipoles a without recourse to the CMB map. In other 
words, our method can be viewed as a form of fossless 
data compression, with all information from the map x 
about the power spectrum coefficient Ci retained in the 
measured coefficients a. 



III. APPLICATION TO WMAP 

Let us now apply our method to the WMAP data [9]. 
After briefly describing the data and foreground masks 
used, we present results for a variety of Galactic cuts. 
To quantify the foreground contribution and select the 
best cut, we then create and analyze simulated CMB and 
foreground maps. This also allows us to optimize other 
practical aspects of our method such as the pixel size. 



A. Data & foreground masks used 

Our analysis is based on the 1 year WMAP data de- 
scribed in [9] with the maps from the five observing 
frequencies combined into the single foreground-cleaned 
map TOH map described in [12] and downloadable at 
http://www.lainbda.gsfc.nasa.gov. The TOH map is 
an all-sky CMB map with the same resolution as the 
highest-frequency WMAP channel (about 12.6' FWHM), 
HEALPix-pixelized at resolution level nside=512. The 
TOH foreground cleaning algorithm assumes that the 
CMB has a blackbody spectrum, but is otherwise com- 
pletely blind, making no assumptions about the CMB 
power spectrum, the foregrounds, WMAP detector noise 
or external templates. 

We explore the same series of increasingly conservative 
foreground masks used in [12] (Figure 1, bottom right). 
Their construction is described in detail in [12], and in- 
volves the following key steps. One first converts the 
maps at the five WMAP frequency bands K, Ka, Q, V 
and W to a common angular resolution and forms four 
difference maps W-V, V-Q, Q-Ka and Ka-K, thereby ob- 
taining maps guaranteed to be free of CMB signal that 
pick up any signals with a non-CMB spectrum. One then 
creates a combined "junk map" (foreground map) by tak- 
ing the largest absolute value of these four maps at each 
pixel. Finally, with appropriate smoothing, one creates 
sky regions based on contour plots of this map. We use 
cuts that are roughly equispaced on a logarithmic scale, 
corresponding to thresholds of 30000, 10000, 3000, 1000, 
300 and 100/xK (Figure 1, right bottom). We label our 
cuts as masks 0, 1, 8, where mask 8 uses pixels with 
less than 100 /xK in the junk map and mask uses the 
entire sky. Masks 1 & 2 merely cut out small blobs in 
the Galactic plane as described in [12]. For masks 0,...,8, 
the total sky percentages removed are 0%, 0.03%, 0.08%, 
0.25%, 0.64%, 2.0%, 6.6%, 24% and 64%, respectively. 




Mask Mask 

FIG. 2: The left and right panels show how the a^^-coefEcients 
(phases) of the quadrupole (left) and the octopole (right) extracted 
from the TOH map change as more sky is masked out. Black 
curves give me the measurements for each mask number, shaded 
bands reflect 1-sigma errors from noise and multipole leakage. For 
example, as explained in detail in Section IIIB, the quadrupole 
error bars (left) do not include quadrupole sample variance the way 
an error bar on a measurement of C2 customarily would, since we 
are interested in the actual a^m-values rather than the underlying 
power spectrum, but they do include sample variance aliased from 
higher multipoles. Based on our simulations, the coefficient most 
susceptible to foreground contamination from the Galactic plane is 
0,20 (top left band), so its mask 6 value is probably closer to the 
truth than its mask value. 



B. Basic results and their interpretation 

For our basic results, we work at HealPix pixel level 
nside = 16 and use the C-matrix corresponding to noise 
and CMB measurements by the WMAP team [9], with 
Ce in equation (3) given by be the best fit ACDM model 
from [1]. As detailed below, our \ow-£ results are quite 
insensitive to these choices. 

Figure 2 shows the result of applying our method to 
this data to measure the components of the quadrupole 
and octopole. The results are plotted as a function of 
Galatic cut, with the thin black lines showing the mea- 
surements aim and the shaded bands illustrating the la 

error range aim i Aaem- These plotted error bars are 

^ ^ 1/2 

AcLijnij — = S^j . The results for more multipoles 

are listed in Table 1. The measurements using the full 

sky (Mask 0) match those reported in [12, 13]. Table 2 

shows the corresponding results for the 3-year WMAP 

data, which was released after the original submission of 

this paper. The results are quite similar — we discuss 
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TABLE I: Measured multipole coefficients in pK. More coefficients are available online. The relation between these real-valued 
spherical harmonic cdcfficituits and the usual complt^x (int^s is given in footiuitc 2. 







Mask 


Mask 5 


Mask 6 


£ 


m 


Fit all 


Fit 0, i,e 


Fit - 3 


Fit - 


- 5 


Fit 0,1, £ 


Fit - 3 


Fit - 5 


2 


-2 


-21.29± 0.20 


-20.48± 1.21 


-20.56± 1.22 


-20.63± 1.22 


-24.07± 2.49 


-24.28± 2.49 


-24.36± 2.53 


2 


-1 


5.94± 0.20 


5.76± 0.23 


5.76± 0.23 


5.76± 0.23 


5.89± 0.40 


5.90± 0.40 


5.93± 0.40 


2 





10.61± 0.22 


13.64± 1.11 


13.69± 1.11 


13.86± 


1.12 


14.04± 2.29 


14.10± 2.30 


14.32± 2.34 


2 


1 


-8.30± 0.20 


-8.21± 0.26 


-8.21± 0.26 


-8.20± 0.26 


-9.14± 0.59 


-9.15± 0.59 


-9.08± 0.60 


2 


2 


-19.39± 0.20 


-20.99± 1.48 


-20.97± 1.48 


-21.16± 1.49 


-17.81± 3.10 


-17.70± 3.11 


-18.53± 3.18 


3 


-3 


40.56± 0.20 


43.17± 1.49 


43.19± 1.50 


43.32± 1.51 


41.84± 3.00 


41.94± 3.00 


42.15± 3.06 


3 


-2 


2.54± 0.20 


2.29± 0.27 


2.29± 0.27 


2.28± 0.27 


2.87± 0.50 


2.86± 0.50 


2.84± 0.50 


3 


-1 


1.00± 0.20 


2.91± 1.03 


2.89± 1.03 


2.78± 1.04 


2.55± 2.10 


2.50± 2.10 


1.86± 2.14 


3 





-6.47± 0.22 


-6.34± 0.29 


-6.34± 0.29 


-6.36± 0.29 


-7.24± 0.65 


-7.23± 0.65 


-7.28± 0.66 


3 


1 


-12.90± 0.20 


-16.45± 1.21 


-16.49± 1.21 


-16.72± 1.22 


-12.05± 2.47 


-12.19± 2.48 


-12.95± 2.55 


3 


2 


30.37± 0.20 


30.76± 0.26 


30.76± 0.26 


30.76± 0.26 


31.53± 0.72 


31.51± 0.72 


31.45± 0.73 


3 


3 


-19.32± 0.20 


-19.02± 1.41 


-18.98± 1.41 


-18.97± 1.42 


-20.76± 2.95 


-20.63± 2.96 


-20.94± 3.04 


4 


-4 


-9.69± 


0.20 


-12.12± 


1.51 




-12.49± 


1.52 


-9.27± 


3.07 




-10.47± 3.10 


4 


-3 


-29.85± 


0.20 


-29.92± 


0.33 




-29.93± 


0.33 


-29.98± 


0.69 




-30.00± 0.69 


4 


-2 


5.24± 


0.20 


4.68± 


1.03 




4.82± 


1.04 


7.65± 


2.07 




7.97± 2.09 


4 


-1 


11.38± 


0.20 


11.66± 


0.28 




11.66± 


0.28 


11.32± 


0.58 




11.26± 0.58 


4 





18.81± 


0.23 


15.95± 


1.10 




15.67± 


1.11 


15.47± 


2.20 




14.87± 2.23 


4 


1 


-8.95± 


0.20 


-9.04± 


0.32 




-9.05± 


0.32 


-7.57± 


0.88 




-7.66± 0.89 


4 


2 


11.75± 


0.20 


13.09± 


1.26 




13.20± 


1.27 


10.98± 


2.52 




11.45± 2.57 


4 


3 


10.23± 


0.20 


9.89± 


0.25 




9.89± 


0.25 


8.42± 


0.74 




8.45± 0.75 


4 


4 


-2.46± 


0.20 


-0.44± 


1.53 




-0.44± 


1.54 


3.64± 


3.10 




3.67± 3.14 


5 


-5 


25.01± 


0.20 


26.42± 


1.61 




26.53± 


1.61 


26.81± 


3.20 




27.09± 3.22 


5 


-4 


11.98± 


0.20 


11.49± 


0.34 




11.48± 


0.34 


11.95± 


0.77 




11.84± 0.77 


5 


-3 


7.70± 


0.20 


5.98± 


1.21 




5.72± 


1.22 


7.05± 


2.34 




6.19± 2.38 


5 


-2 


2.46± 


0.20 


2.78± 


0.31 




2.78± 


0.31 


2.00± 


0.66 




2.04± 0.66 


5 


-1 


5.31± 


0.21 


3.46± 


1.01 




3.49± 


1.01 


4.31± 


2.01 




4.46± 2.02 


5 





15.61± 


0.23 


15.34± 


0.34 




15.35± 


0.34 


16.46± 


0.83 




16.44± 0.83 


5 


1 


30.49± 


0.21 


34.32± 


1.18 




34.40± 


1.18 


31.00± 


2.30 




31.23± 2.35 


5 


2 


-11.21± 


0.20 


-11.76± 


0.32 




-11.78± 


0.32 


-12.68± 


0.95 




-12.73± 0.96 


5 


3 


26.11± 


0.20 


25.84± 


1.15 




25.82± 


1.16 


26.67± 


2.28 




26.81± 2.33 


5 


4 


-6.84± 


0.20 


-7.14± 


0.28 




-7.15± 


0.28 


-5.49± 


0.83 




-5.55± 0.83 


5 


5 


14.00± 


0.20 


12.74± 


1.53 




12.76± 


1.54 


14.22± 


3.10 




14.29± 3.12 



the differences below in Section IV and Section V.^ The 
quadrupole moments were computed by taking a to in- 
clude all multipole coefRcients with £ < 2. It is important 
to always include the monopole {£ = 0) and dipole = 1) 
so that our method makes the measurements that we are 
interested in (in this case £ = 2) completely independent 
of these totally unknown quantities (the maps released 
by the WMAP team have already been approximately 



^ We foreground-cleaned the 3 year WMAP data with the exact 
same TOH algorithm [12] as was used for the 1-year data, and 
the interested reader can download the cleaned map from 
http : //space .mit . edu/home/tegmark/wmap . html. 



purged of a 3K monopole corresponding to the mean 
CMB temperature and a 6mK dipole corresponding to 
Earth's motion). Similarly, the octopole measurements 
were computed by taking a to include all multipole coef- 
ficients with £ = 3 and £ <1 (We will use the alternative 
£ < 3 option in Section IV where we study quadrupole- 
octopole alignment and want uncorrelated measurements 
of the two.) 

Note that the error bars Aaim on, say, the quadrupole 
coefRcients shown in Figure 2 do not include quadrupole 
sample variance the way an error bar on a measurement 
of C2 customarily would. This is why they are seen to 
be so tiny for Mask 0, reflecting only the detector noise 
contribution. The error bar Aa£,„ simply reflects the 
uncertainty in our measurement of the coefficient aim 
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TABLE II: Same as previous table, but for the 3-year WMAP data. 







Mask 


Mask 5 


Mask 6 


I 


m 


Fit all 


Fit 0, l,i 


Fit - 3 


Fit 0- 


- 5 


Fit 0, 1,1 


Fit - 3 


Fit - 5 


2 


-2 


-24.29± 0.20 


-23.06± 1.21 


-23.14± 1.22 


-23.20± 1.22 


-25.88± 2.49 


-26.09± 2.49 


-26.20± 2.53 


2 


-1 


6.83± 0.20 


6.70± 0.23 


6.69± 0.23 


6.69± 0.23 


6.65± 0.40 


6.66± 0.40 


6.69± 0.40 


2 





3.22± 0.22 


3.88± 1.11 


3.94± 1.11 


4.11± 1.12 


3.79± 2.29 


3.86± 2.30 


4.14± 2.34 


2 


1 


0.37± 0.20 


0.51± 0.26 


0.51± 0.26 


0.52± 0.26 


-0.62± 0.59 


-0.62± 0.59 


-0.54± 0.60 


2 


2 


-21.06± 0.20 


-19.54± 1.48 


-19.52± 1.48 


-19.72± 1.49 


-15.31± 3.10 


-15.18± 3.11 


-16.11± 3.18 


3 


-3 


45.13± 0.20 


47.03± 1.49 


47.06± 1.50 


47.24± 


1.51 


44.50± 3.00 


44.59± 3.00 


44.97± 3.06 


3 


-2 


2.36± 0.20 


2.11± 0.27 


2.11± 0.27 


2.10± 0.27 


2.62± 0.50 


2.62± 0.50 


2.59± 0.50 


3 


-1 


8.44± 0.20 


10.58± 1.03 


10.55± 1.03 


10.46± 


1.04 


8.72± 2.10 


8.66± 2.10 


8.07± 2.14 


3 





-8.45± 0.22 


-8.27± 0.29 


-8.27± 0.29 


-8.28± 0.29 


-9.52± 0.65 


-9.50± 0.65 


-9.54± 0.66 


3 


1 


-15.88± 0.20 


-15.97± 1.21 


-16.00± 1.21 


-16.24± 1.22 


-11.94± 2.47 


-12.03± 2.48 


-12.86± 2.55 


3 


2 


32.65± 0.20 


32.96± 0.26 


32.96± 0.26 


32.96± 0.26 


33.71± 0.72 


33.70± 0.72 


33.61± 0.73 


3 


3 


-16.76± 0.20 


-18.59± 1.41 


-18.55± 1.41 


-18.52± 1.42 


-20.71± 2.95 


-20.60± 2.96 


-20.79± 3.04 


4 


-4 


-12.24± 


0.20 


-14.39± 


1.51 




-14.81± 


1.52 


-11.83± 


3.07 




-13.10± 3.10 


4 


-3 


-27.82± 


0.20 


-27.77± 


0.33 




-27.79± 


0.33 


-27.60± 


0.69 




-27 62± 69 


4 


-2 


5.74± 


0.20 


4.79± 


1.03 




4.92± 


1.04 


7.10± 


2.07 




7.43± 2.09 


4 


-1 


9.79± 


0.20 


9.99± 


0.28 




lO.OOi 


0.28 


9.89± 


0.58 




9.84± 0.58 


4 





17.25± 


0.23 


16.71± 


1.10 




16.43± 


1.11 


16.62± 


2.20 




16.02± 2.23 


4 


1 


-14.43± 


0.20 


-14.63± 


0.32 




-14.64± 


0.32 


-12.94± 


0.88 




-13.02± 0.89 


4 


2 


16.33± 


0.20 


15.01± 


1.26 




15.12± 


1.27 


12.45± 


2.52 




12.89± 2.57 


4 


3 


12.17± 


0.20 


11.91± 


0.25 




11.90± 


0.25 


9.89± 


0.74 




9.92± 0.75 


4 


4 


-5.72± 


0.20 


-3.22± 


1.53 




-3.23± 


1.54 


1.14± 


3.10 




1.19± 3.14 


5 


-5 


27.10± 


0.20 


28.39± 


1.61 




28.52± 


1.61 


28.69± 


3.20 




29.01± 3.22 


5 


-4 


11.39± 


0.20 


10.89± 


0.34 




10.88± 


0.34 


11.41± 


0.77 




11 28± 77 


5 


-3 


3.89± 


0.20 


2.69± 


1.21 




2.40± 


1.22 


4.57± 


2.34 




3.63± 2.38 


5 


-2 


2.56± 


0.20 


2.89± 


0.31 




2.90± 


0.31 


2.21± 


0.66 




2.26± 0.66 


5 


-1 


0.90± 


0.21 


-1.15± 


1.01 




-1.14± 


1.01 


0.98± 


2.01 




1.09± 2.02 


5 





15.78± 


0.23 


15.45± 


0.34 




15.45± 


0.34 


16.98± 


0.83 




16.94± 0.83 


5 


1 


34.10± 


0.21 


34.50± 


1.18 




34.58± 


1.18 


31.71± 


2.30 




31.94± 2.35 


5 


2 


-9.74± 


0.20 


-10.19± 


0.32 




-10.21± 


0.32 


-11.06± 


0.95 




-ll.lOi 0.96 


5 


3 


24.12± 


0.20 


25.59± 


1.15 




25.55± 


1.16 


26.41± 


2.28 




26.53± 2.33 


5 


4 


-5.72± 


0.20 


-5.95± 


0.28 




-5.97± 


0.28 


-3.74± 


0.83 




-3.83± 0.83 


5 


5 


13.47± 


0.20 


12.37± 


1.53 




12.41± 


1.54 


13.25± 


3.10 




13.37± 3.12 



on the sky seen from our particular vantage point in 
space. These coefficients vary randomly with a variance 
C2 between widely separated Hubble volumes, and this 
quadrupolc sample variance therefore enters only in the 
next data analysis step where our five locally measured 
quadrupole coefficients are used to make inferences about 
the value of C2. 

As mentioned in Section II, however, the error bars 
Aa^TO do include sample variance from the other multi- 
poles that were excluded from the a- vector, i.e., i > 3 
for our quadrupolc example. This is why the error bars 
in Figure 2 are seen to flare up dramatically towards the 
right: as large portions of the sky get masked out, the 
severely broken orthogonality between the spherical har- 
monics makes higher multipoles contribute to quadrupole 



estimators. Since the values of these higher multipole 
coefficients are unknown to us and are not measured by 
the method in this case, their contribution is counted as 
noise which can be quantified since their variance C{ 
is known. Since our method tries to minimize the error 
bars, it makes the best tradeoffs it can to minimize the 
net effect of such leakage from non-included mTiltipolcs. 
For multipoles included in a, in contrast, our method 
forces the leakage to be exactly zero, at the cost of larger 
error bars and more leakage from non-included multi- 
poles. 

The disk-like geometry of our Galaxy provides some 
intuition for the behavior of the error bars. Because the 

Galaxy cuts are approximately symmetric under reflec- 
tion (parity-even), different values couple mainly if they 
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arc separated by an even number. To the extent that the 
Galaxy cut is azimuthaUy symmetric (a crude approx- 
imation at best), different m- values do not couple. For 
example, the quadrupole estimator 021 therefore picks up 
"noise" mainly from 041 , aei , etc. Moreover, since spheri- 
cal harmonics with high \m\ live mainly near the Galactic 
equator and rapidly approach zero towards the Galactic 
poles, they are the ones that suffer most as the cut is in- 
creased and therefore have the largest error bars in Table 
1. 

The above-mentioned sample variance considerations 

also shed light on how the method is affected by changes 
in the C-matrix, i.e., on how robust our method is to as- 
sumptions about unwanted noise and CMB signals. Since 
the method is required to faithfully (without bias) recover 
the multipoles a whatever they are, one would intuitively 
expect that W and therefore a are independent of con- 
tributions from signal and noise in these multipoles to 
the C-matrix. With some matrix algebra, one can prove 
that that this is indeed the case, specifically that adding 
to C a contribution of the form YS5Y* for some ma- 
trix S5 leaves W and a unchanged and simply increases 
S by S5. For example, if we use our method to mea- 
sure the five quadrupole coefhcients, the result a will be 
independent of our assumptions about the true value of 
C2. 



C. Simulations with foregrounds 

The key remaining issue is how to quantify the contri- 
bution of residual foreground contamination and how to 
optimize the Galaxy cut to minimize this contribution. 
We will now address this issue with simulations. 

The residual foreground contamination present after 
foreground cleaning can be quite important. All five lin- 
early cleaned WMAP-maps used in the literature have 
their problems. The WMAP team ILC map [9] comes 
with the disclaimer that it should not be used for scien- 
tific analysis — needless to say, the science questions are 
so interesting that this has not stopped large numbers of 
groups from using it anyway, particularly when all-sky 
coverage was helpful. The TOH cleaned WMAP map 
[12] has the advantage of having lower total noise and 
residual foregrounds, but like the ILC map, it lacks er- 
ror bars accurately quantifying these residuals, as do all 
other foreground-cleaned maps published to date. The 
SMICA map [108] is generated with a method very si- 
milar that used for the TOH map [12, 107], and the re- 
sult is indeed encouragingly similar [108]. The WI-FIT 
map of [109] is also closely related (since this method 
cleans with linear combinations of linear combinations 
(pairwise differences) of maps, it too is a linear and 
WMAP-internal technique just like the ILC, TOH and 
SMICA methods — the only fundamental difference be- 
tween the four methods is that the approximate variance 
minimizes minimization is performed in pixel space for 
ILC, in (masked) spherical harmonic space for TOH and 



SMICA and in wavelet space for WI-FIT (and optionally 
SMICA) . Finally, using external template maps alone [9] 
suffers from both dubious extrapolation and from the 
problem that they {e.g., the Haslam- based synchrotron 
map) have striping or other systematic errors at a sub- 
stantially larger level than the WMAP maps, and that 
these errors will propagate into the cleaned map, perhaps 
making it less clean that it was to start with in the most 
systematics-affected modes — template-cleaned residual 
WMAP foregrounds have never been accurately quanti- 
fied in the literature. 

A useful recent approach has been to quantify ILC er- 
rors by reapplying the ILC cleaning algorithm to simu- 
lated sky maps where the correct answer is known and the 
residuals can thus be computed [110]. This has confirmed 
that the residual foreground contribution is substantial. 
We will adopt a similar procedure for the TOH map, but 
with an important difference that makes our approach 
more conservative. 



1. Mock foregrounds 

It is absolutely crucial to use a foreground model that 
does not artificially mimic the cleaning algorithm. For 
example, as acknowledged in Eriksen et al. (2005) [110], 
their simulations underestimate the residual foregrounds 
by assuming that the frequency dependence of each phys- 
ical component is the same in all directions in the sky — 
this assumption of perfect frequency coherence [111] is 
empirically known to be incorrect for synchrotron radi- 
ation (see, e.g., [112] for a recent treatment), and al- 
lows the simulated foregrounds to be perfectly removed 
by the cleaning method. To avoid reaching overly op- 
timistic conclusions, we therefore make simulated fore- 
ground templates at the five WMAP frequencies in a 
different way: straight from data by subtracting the 
TOH CMB maps from the five WMAP frequency chan- 
nels. The resulting five foreground maps that we use are 
shown in in Figure 1 of [12]. This procedure is prob- 
ably overly conservative because some foregrounds and 
noise get double-counted: for instance, our K-band "fore- 
ground" map includes detector noise from all five WMAP 
frequencies because it was constructed by subtracting the 
TOH map (with its noise) from the K-band observations. 



2. Mock CMB 

We generate a simulated CMB-only sky map (Figure 1, 
top right) by using the HEALPix software and the best fit 
WMAP power spectrum from [1], smoothed to have the 
same angular resolution as the TOH map (corresponding 
to the beam of the W-band, the WMAP channel with the 
highest resolution). As input to the cleaning algorithm, 
we also smoothed it to the four lower resolutions corre- 
sponding to the K-, Ka-, Q- and W-band beams [113]. 
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3. Mock detector noise 

We construct simulated WMAP noise maps at the five 
channels by assuming that the noise is Gaussian and un- 
correlated both between different pixels and between dif- 
ferent channels, which is a good approximation accord- 
ing to [9]. We set the noise amplitude in each map so 
that the rms fluctuation levels match those estimated 
from the WMAP team (Table 5 in [114]). We con- 
firmed that this noise normalization matched both the 
rms in various same-frequency difference maps such as 
W = {Wl -W2 + W3- W4)/4 and the total power 
spectrum normalization at very high £ where CMB and 
foregrounds are negligible. Specifically, we used rms 
noise values N|/^ of 150/^K, 132^K, 122^K, U9^iK and 
179/iK for K, Ka, Q, V and W, respectively. This cor- 
responds to an ^-independent noise power spectrum of 
order Cg ^ 0.1/iK^ in each channel, i.e., orders of mag- 
nitude below the large-scale CMB signal ([12] estimates 
C2 200AiK2 and C2 ~ SOO^iK^), so the details of our 
noise modeling have no impact whatsoever on our re- 
sults — rather, our uncertainties are dominated by sam- 
ple variance and residual foregrounds. We tried both 
crude noise maps with equal variance in all pixels and 
more accurate ones with the variance modulated in in- 
verse proportion to the number of WMAP observations 
of each pixel. As expected from the fact that noise is 
negligible at the low ^-values of interest here, these two 
approaches give for all practical purposes indistinguish- 
able results — we adopt the former for simplicity. 

4. Residuals 




2 4 6 8 

Mask 



FIG. 3: The left and right panels show how the measurement 
errors Aaf,„ = atrn ~ "^f^" quadrupole (left) and 

octopole (right) extracted from a simulated map change with in- 
creasing Galactic cut. The three curves in each panel correspond 
to errors from CMB multipole leakage alone (black) CMB plus de- 
tector noise (red) and CMB plus noise plus residual foregrounds 
(blue) The shaded/green areas show the la error bars correspond- 
ing to leakage and noise. Note that these foreground-related er- 
rors are typically an order of magnitude smaller than the expected 

1/2 

CMB signal (the ACDM concordance model gives 
1/2 

and Cg = 20/xK), and that they are even more subdominant for 
power spectrum estimation where they add only in quadrature. 



For each of the five WMAP channels, we sum the mock 
CMB, noise and foreground maps described above. We 
then apply the TOH foreground cleaning method to these 
five input maps exactly as described in [12] to obtain 
a single foreground-cleaned map. Finally, we subtract 
the "true" simulated CMB sky map from the cleaned 
map, obtaining the residual map shown in Figure 1 (bot- 
tom left). This figure visually illustrates that the TOH 
method provides an unbiased recovery of the CMB, since 
the residual map shows no patterns whatsoever in com- 
mon with the input CMB map (upper right). Compar- 
ing the residual map with the V-band map (top left) also 
illustrates that the TOH method reduces the levels of 
both foregrounds and detector noise. Nonetheless, resid- 
ual foregrounds are clearly visible, particularly near the 
Galactic plane in the parts that the "junk map" (bottom 
right) suggests are highly contaminated. 

If we knew that this was the actual residual in our TOH 
map, we could of course eliminate the foreground prob- 
lem completely by subtract it out. This is unfortunately 
not the case. To be conservative, will make no attempt 
to further clean our input map, and will merely use the 
residual map to quantify the errors in our multipole mea- 
surements. 



What does the residual map depend on? By repeating 
the entire procedure with independent CMB simulations, 
we found that the residual map was rather independent 
of what CMB map was used. This is expected, since the 
TOH method is completely CMB-independent except for 
the small effect that CMB sample variance can have on 
the computation of optimal weights for the cleaning pro- 
cess. Repeating the entire procedure with independent 
noise simulations alters the residual map mainly on the 
very smallest angular scales. The large-scale patterns 
visible in Figure 1 (bottom left) are therefore determined 
mainly by properties of the foregrounds themselves. For 
example, foregrounds whose spectral index varies notably 
across the sky cannot be perfectly subtracted out by the 
TOH method. Since these large-scale residual patterns 
are all that matter for our present focus on low there 
is no point in making more than one simulation of our 
type, i.e., it is pointless to make an ensemble of simula- 
tions with different CMB and noise realizations. 
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D. Optimizing the Galaxy cut 

To quantify the effect of residual foregrounds on our 
multipole measurements, we repeat tlie analysis from 
Section III using our simulated maps. Specifically, we 
measure the Ofm-coefficients from the simulated CMB- 
only map, from a map with CMB and noise (made with- 
out introducing any foregrounds before the cleaning pro- 
cess) and from the cleaned map including CMB, noise 
and foregrounds. The results are shown in Figure 3 with 
the correct (simulated) a^m-values subtracted off. 

The fact that the black and red curves are almost iden- 
tical confirms the above-mentioned claim that WMAP 
detector noise makes essentially no difference on these 
large angular scales. The fact that the red curves 
agree well with the shaded regions confirms that our 
method and our software implementing it are working 
as they should, i.e., that the multipoles are faithfully 
recovered with accuracy consistent with the predicted 
noise-l- leakage error bars Aa^„j . 

The key new information in Figure 3 is contained in 
the blue curves, which include errors caused by residual 
foregrounds. These curves reveal two interesting facts: 

1. Residual foreground contamination near the Galac- 
tic plane afflicts the a2o-component at a level 
around ten /iK. 

2. Residual foreground contamination at higher lati- 
tudes appears to afflict most coefficients at the level 
of a few /iK. 

The fact that 020 is special follows from the well-known 
fact that this is the component most similar in shape 
to the Galactic emission, as pointed out already in the 
COBE/DMR analysis [8]. Specifically, in the very crude 
approximation that the Galactic plane has both parity 
and azimuthal symmetry, it would contaminate only the 
multiple coefficients 020, 040, ^60: e^c. Figure 3 shows 
that most of this residual contamination comes from the 
innermost parts of the plane, so that cutting out merely 
the dirtiest 7% of the sky (Mask 6) cuts the foreground 
contribution to 020 from about 9/iK to 3/iK. 

In contrast, and in good agreement with these symme- 
try considerations, none of the other quadrupole or oc- 
topole coefficients show evidence of contamination from 
the inner Galactic plane: the foreground contributions 
shown by the blue curves do not tend to shift closer to 
the zero as the cut is increased up to Mask 6. Beyond 
this, the leakage errors are seen to exceed the residual 
foreground errors. This small residual foreground conta- 
mination is seen to typically be 1 — 3/iK, of similar mag- 
nitude to the a2o-contamination that remains for Mask 
6. 

So based on these results, what is the optimal galaxy 
cut when measuring multipole coefficients? We have 
seen that too small a Galaxy cut leads to unneccessar- 
ily high foreground contamination, whereas too large a 
cut leads to unneccessarily high leakage from other mul- 
tipoles. The natural optimum is seen to be Mask 5 or 




Mask 

FIG. 4: Effects of pixelization. The panels show how the measure- 
ment errors Aa^^ = — a^l^" in /iK extracted from a simulated 
map change with increasing Galactic cut. the black, blue and red 
curves correspond to Healpix pixelization at angular resolution level 
•aside =4, 8 and 16, respectively. The shaded/green bands show 
the Icr errors from detector noise and leakage. 

6, where the typical errors Aa^m are of order a /iK from 
leakage and a few /iK from foregrounds, and we have 
therefore listed results for both in Table 1. For the more 
aggressive Mask 7, sample variance strongly dominates, 
and for the less aggressive Mask 4, the 020 foreground 
contamination is noticeably larger. If the reader wishes 
to perform an analysis that it no way involves leakage 
from the black sheep 020, another natural choice is Mask 
0. 



E. Optimizing the pixel size 

This subsection discusses practical details useful for 
readers interested in applying this method. Other readers 
may wish to skip straight to Section IV. 

Applying equations (4) and (5) using the native 
WMAP pixel resolution {nside = 512) is unfeasible, since 
it would involve the inversion of an A'^ x A'^-dimensional 
C-matrix with A^ of order 3 million. Since we are only 
interested in the lowest multipoles, there are fortunately 
two simple ways to overcome this problem. The first 
approach involves using matrix identities to transform 
the problem into one involving smaller matrices (as for 
the so-called Woodbudy formula). The second approach, 
which we used above, is to simply smooth the map onto 
a coarser pixelization and correct measurements for this. 

Figure 4 shows the errors on the recovered quadrupole 



10 




FIG. 5: Change in the direction of the quadrupole and octopole "axis of evil" directions as a function of the mask number. The two 
crosses correspond to the true quadrupole and octopole axes in our CMB simulation, chosen to match the WMAP measured values. The 
departures of the shaded dots from the corresponding crosses reflect the effect of noise, multipole leakage and foreground contamination. 



coefficients from our simulated map (witfi neither detec- 
tor noise nor foregrounds) for pixelizations with nside—A, 
8 and 16, showing that as long as we work at resolution 
16, pixelization errors are negligibly small compared to 
noise and sample variance for the masks of interest. 

This pixelization corresponds to averaging the under- 
lying sky map across the area of each pixel. If the map 
were subjected to Gaussian smoothing, then the multi- 
poles aim would be suppressed by the well-known factor 
^eH{e+i)/2^ where 9 =FWHM/\/ihr2 and FWHM is the 
full-width half-maximum of the Gaussian smoothing ker- 
nel. For our adopted resolution level 16, the pixels are 
approximately squares of side 3.7°. A very crude estimate 
therefore suggests that the corresponding multipole sup- 
pression is of order e~(€/40)^/2, i.e., a 0.3% quadrupole 
suppression and a 2% suppression for £ = 5. Since this 
smoothing effect is negligible compared to the measure- 
ment errors caused my leakage and foregrounds, we have 
not corrected for it in Table 1. 

In order to quantify the errors due to pixelization, we 
calculate the five components of the quadrupole for dif- 
ferent nsides of 4,8 and 16, and also studied how these 
values change as we increase the mask number (or in- 
crease the Galactic cut). 

Although resolution 8 already gives us a good estimate 
of the lower multipoles, we decided (to be conservative) 
to use nside=16 in all calculations done below. 



IV. IMPLICATIONS FOR THE 
QUADRUPOLE-OCTOPOLE ALIGNMENT 

As a simple example of an application of our results, 
we use our new multipole measurements to revisit the sig- 
nificance of the WMAP quadrupole-octopole alignment. 
The two crosses in Figure 5 show the "axis of evil" from 
[12, 13], extracted from a simulated CMB map where 
we inserted the exact same multipole and quadrupole 
coefficients as measured from WMAP. We then added 
noise and foregrounds as above and extracted the "axis 
of evil" for our galactic cut. Confirming the conclusion of 
[28], the octopole is seen to be quite robost, whereas the 
quadrupole moves around somewhat more (it is clearly 
more fragile due to its intrinsically lower amplitude). As 
seen in Figure 6, this causes the apparent measured align- 
ment be somewhat less significant than the true one, but 
makes no dramatic difference. Similarly, we find that re- 
placing Mask by Mask 6 (the joint quadrupole/octopole 
fit in Table 1) degrades the alignment significance only 
slightly, from a one-in-sixty fluke to a one-in-forty fluke. 



V. CONCLUSIONS 

We have presented a minimum-variance method for 
measuring the CMB multipole coefficients aim given 
anisotropic noise, incomplete sky coverage and fore- 
ground contamination. Our method constitutes loss- 
less data compression in the sense that the widely used 
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FIG. 6: The top panel shows the separation angle (in degrees) be- 
tween the quadrupole and octopole vectors and how they change as 
a the mask number increases. The lower panel shows the probabil- 
ity for this alignment to happen. In both figures the red, blue 
and green lines represent the coefficient values for a map with 
CMB only, CMB plus noise, and CMB plus noise plus residual 
foreground, respectively. 

quadratic estimators of the power spectrum Ci can be 
computed directly from our Ofm-estimators. We illus- 



trated the method by applying it to the WMAP data. 
As the Galactic cut is increased, the error bars Aaim on 
low multipoles go from being dominated by foregrounds 
to being dominated by sample variance from other multi- 
poles, with the intervening minimum defining the optimal 
cut. 

Because WMAP detector noise is negligible on these 
large angular scales, one would not expect the improved 
sensitivity of second WMAP data release to significantly 
modify our results per se. Comparing Table 1 with Ta- 
ble 2 shows that the multipoles that change significantly 
are those most sensitive to foregrounds, notable 020, sug- 
gesting the interpretation that the lower noise levels im- 
proved the foreground cleaning. The corresponding mask 
6 quadrupole and octopole "axes of evil" change from 
{I, h) = (251.5°, 52.1°) and (238.4°, 63.0°) for WMAPl to 
{l,h) = (222.4°, 75.5°) and (236.4°, 64.9°) for WMAP3, 
respectively. In other words, the octopole axis is essen- 
tially unchanged whereas the quadrupole axis moves to 
the other side of the octopole axis, getting still closer 
(from 13.0° off to 11.5°), requiring a l-in-50 fluke - see 
Figure 7. Moving the quadrupole axis in a random di- 
rection would typically degrade the alignment, so this is 
consistent with the hypothesis that the true alignment 
is still better and that this is partially masked by fore- 
grounds. 

Improved sensitivity from WMAP and Planck and also 
new cleaning techniques will hopefully allow better re- 
moval and quantification of residual foreground contami- 
nation, which would reduce the main source of uncer- 
tainty in our measurements and shed further light on the 
low-^ CMB puzzles. 
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FIG. 7: Changes in the direction of the preferred quadrupole and octopole axes as a function of mask number for the WMAP 1-year 
(left) and 3-year (right) data. In both panels, the quadrupole and octopole results are represented by balls and stars, respectively, and the 
colors black, cyan, red, blue, magenta, green and yellow correspond to the results for masks from to 6, respectively. The quadrupole axis 
is seen to be more robust in the 3-year data, suggesting that the lower noise level may have improved the performance of the foreground 
cleaning process, reducing the level of residual low-latitude foreground contamination. 
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